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1  NORSAR  records  of  P  and  P  cods  from  STRAIT  and  POOL  5 

showing  time  windows  selected  for  Initial  P  and  P  coda. 

Note  that  the  energy  in  the  initial  P  for  the  two 
explosions  with  nearly  equal  yields  is  considerably 
different  whereas  the  P  coda  levels  do  not  appear  to  be 
much  different. 

2a  Results  from  STRAIT  recorded  at  NORSAR.  Displacement  6 

spectra,  not  corrected  for  Instrumental  response,  of  the 
observed  Initial  P  signal  (S)  denoted  by  the  symbol  +, 
noise  (N),  symbol  0,  for  a  window  of  equal  length,  and 
S-N,  symbol  *. 

2b  Results  from  STRAIT  recorded  at  NORSAR.  Displacement  7 

spectra  of  P  coda,  similar  to  those  in  Figure  2a. 

2c  Results  from  STRAIT  recorded  at  NORSAR.  Spectral  ratio  R 

P/P-coda  (corrected  for  noise),  points  for  which  the 
signal  to  noise  power  ratio  Is  leas  than  2  are  not 
plotted. 

3a  Results  from  MAST  recorded  at  NORSAR,  similar  to  those  10 

for  STRAIT  In  Figure  2a.  Spectral  nulls  In  Initial  P  at 
about  1  and  2  Hz  are  distinct  and  can  also  be  seen  in 
the  spectral  ratio  plot  of  Plgure  3c. 

3b  Results  from  MAST  recorded  at  NORSAR,  similar  to  those  11 

for  STRAIT  in  Figure  2b. 

3c  Results  from  MAST  recorded  at  NORSAR,  similar  to  those  12 

for  STRAIT  In  Figure  2b. 

4  Mean  P/P-coda  spectral  slope  for  the  frequency  range  0.5  13 

to  3.0  Hz  versus  average  ahot-polnt  to  surface  or  medium 
velocity  for  14  explosions.  The  linear  trend  (dashed 

line)  shows  correlation  coefficient  of  0.R96. 

5  Spectral  particle  velocity  (arbitrary  units  and  integrated  lfi 

over  the  frequency  range  0.5  to  3.0  Hz)  plotted  against  log 

yield  for  the  Initial  P  from  the  fourteen  NTS  explosions. 

Lines  of  slope  of  0.9  are  drawn  through  data  points  for  the 
four  Yucca  Flats  events  (upper  line)  and  the  ten  Pahute  Mesa 
events  (lower  line).  For  actual  data  points,  see  Appendix 
A. 

6  Same  as  Figure  5  for  P  coda.  A  single  line  of  slope  0.9  19 

is  drawn  through  data  points  for  all  fourteen  events. 

For  actual  data  points,  see  Appendix  A. 
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7  Initial  P  Integrated,  spectral  displacements  with  t*  and  20 

granite  source  function  corrections  (arbitrary  units)  versus 

log  yield  for  the  fourteen  NTS  explosions.  Lines  of  slope 
1.0  are  drawn  separately  for  the  four  Yucca  Flats  events 
(upper  line)  and  the  ten  Pahute  Mesa  events  (lower  line). 

If  a  single  line  of  slope  unity  were  drawn  through  all  data 
points,  the  standard  deviation  of  the  intercept  would  be 
0.14  magnitude  unit.  For  actual  data  points,  see  Appendix 
A. 

8  P  coda  integrated  spectral  displacement  with  t*  and  granite  21 

source  function  corrections  (arbitrary  units)  plotted  as  a 
function  of  log  yield.  A  line  of  slope  1.0  through  all  data 

points  has  an  associated  standard  deviation  of  0.06 
magnitude  unit.  For  actual  data  points,  see  Appendix  A. 

9  P  coda  Integrated  spectral  displacement  with  t*  and  tuff  22 

source  function  corrections  (arbitrary  units)  plotted  versus 

log  yield.  A  line  of  slope  1.0  through  11  data  points 
(excluding  the  one  with  the  smallest  yield)  gives  a  standard 
deviation  of  0.07  magnitude  unit.  For  actual  data  points, 
see  Appendix  A. 

10  Similar  to  Figure  9  except  that  the  tuff  source  function  21 

correction  Is  based  on  approximate  yields  derived  from  m^ 
values.  A  line  of  slope  1.0  through  the  thirteen  data 

points  Indicates  a  standard  deviation  of  0.09  magnitude 
unit.  For  actual  data  points,  see  Appendix  A. 

11  Initial  P  integrated  spectral  displacement  with  t*  and  28 

granite  source  function  corrections  (arbitrary  units) 

plotted  as  a  function  of  m.  for  14  Soviet  explosions.  Mean 
slope  Is  0.75  +  0.05  and  tne  correlation  coefficient  is 
0.97.  ~ 

12  P  coda  Integrated  spectral  displacement  with  t*  and  granite  29 

source  function  corrections  (arbitrary  units)  versus  hl  for 

23  Soviet  explosions.  Mean  slope  is  0.70  +  0.03  and  tne 
correlation  coefficient  is  0.98. 

13  P  coda  spectral  displacement  (arbitrary  units)  plotted  31 

versus  yield  values  obtained  from  Lg  (Nuttli,  1983)  on  a 

log-log  scale.  The  correlation  coefficient  is  0.82. 
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INTRODUCTION 


Scattering  of  selaalc  waves  due  to  departure  froa  plane  stratified 
aedia  appears  to  play  a  vital  role  in  the  observed  ground  notion  froa 
nuclear  explosions.  A  striking  exaaple  is  the  transverse  component  Lg 
which  la  generally  larger  than  the  vertical  or  the  radial  component 
(Gupta  and  Blandford,  1983).  Transverse  Lg  froa  explosions  at  the 
Nevada  Test  Site  (NTS)  has  been  found  to  be  alao  a  better  estimator  of 
yield  than  direct  P  (von  Seggern  and  Alexander,  1984).  At  teleselaaic 
distances,  algnif icantly  greater  stability  of  P  coda  (amplitudes  and 
spectra)  as  compared  to  direct  P  has  been  pointed  out  and  established  in 
nuaerous  recent  studies  (see,  e.g.  the  review  articles  by  Akl,  1982  and 
Bache,  1982).  For  a  Halted  test  area,  P  coda  should  therefore  provide 
a  aore  precise  aeaaureaent  of  yield  than  the  initial  P.  There  are 
several  nechanlsae  that  aay  operate  to  aake  coda  amplitudes  aore  precise 
than  the  primary  signal.  P  coda  consists  of  scattered  energy  where  most 
of  the  scattering  takes  place  near  the  source  and  the  receiver. 
Near-source  scattering  is  perhaps  mainly  due  to  the  scattering  of 
source-generated  surface  waves  to  P  body  waves  by  topographic  features 
and  other  iapedance  irregularities  within  several  klloaeters  of  the 
source  (Greenfield,  1971).  Contribution  to  teleselaaic  P  froa 
near-eource  scattering  occurs  for  signals  emerging  over  a  large  portion 
of  the  focal  sphere  so  that  the  effects  of  reflection  froa  the  free 
surface  (i.e.,  cancellation  by  pP)  will  be  suppressed.  Among  others, 
Blandford  and  Shuaway  (1982)  have  remarked  on  the  large  variability  of 
teleselaaic  P  due  to  the  effect  of  pP  so  thst  the  use  of  coda  aay  well 
reaove  this  problem.  The  uniformity  of  coda  amplitudes  over  the 
eleaents  of  an  array  for  a  single  event  suggests  that  coda  is 
insensitive  to  the  near-receiver  focusing  and  defocusing  effects  that 
can  drastically  Influence  the  primary  arrival.  This  probably  happens 
because  (1)  the  energy  scattered  near  the  source  and  arriving  directly 
at  the  receiver  travels  not  along  the  path  followed  by  the  direct  P  but 
along  nuaerous  adjacent  paths  so  that  the  coda  aay  represent  an  average 
over  the  amplitudes  generated  within  a  relatively  large  source  region; 
and  (2)  direct  signals  are  received  In  a  relatively  large  region 
surrounding  the  receiver  and  then  scattered  to  the  receiver,  thus 
receiver  focusing  is  averaged  out. 


Local  or  near-recelver  scattering  should  Influence  lnlttal  P  ae 
well  aa  the  later  arrivals.  In  fact,  topography  and  other 
lnhoaogeneltles  near  the  receiver  should,  on  the  hanls  of  reciprocity, 
be  aa  significant  scatterers  as  they  are  when  near  the  source.  Oupta 
and  Blandford  (1981)  analyzed  the  spectral  composition  of  the  three 
orthogonal  components  of  ground  motion  in  the  Initial  P  waves  from 
SALMON  and  concluded  that  local  scattering  Is  significant,  especially  at 
high  frequencies.  One  may,  however,  expect  near-source  scattering  to 
dominate  over  local  scattering  In  some  cases  and  local  scattering  to 
greatly  exceed  near-source  scattering  in  certain  other  situations.  Tn 
other  words,  only  certain  "clean**  sites  may  be  suitable  for 
Investigating  near-source  scattering.  In  this  study,  observations  from 
the  larg<i  aperture  seismic  array,  NORSAR,  located  on  Permian  and  older 
hard  (mostly  gneisses  and  granite)  rocks  are  used.  Large  subsurface 
impedance  contrasts  should  not  exist  at  this  locality  so  that 
near-recelver  scattering  should  be  low.  Most  of  our  analyses  of  P  and  P 
coda  Is  carried  out  on  records  of  explosions  at  NTS  for  which  the  yield 
and  shot  medium  properties  were  known.  Records  of  Soviet  explosions  at 
the  Shagan  site  In  the  Semlpalatlnak  region  are  also  examined  and  their 
yields  estimated  on  the  basis  of  the  NTS  experience. 


COMPARISON  OP  P  AND  P  CODA  SPP.CTRA 


Short-period ,  vertical  component  NORSAR  records  of  fourteen  (14) 
underground  nuclear  explosions  (Table  l)  were  first  selected  for 
spectral  analysis  of  initial  P  and  the  associated  codAS.  The  selection 
was  aade  so  as  to  observe  the  following  criteria:  (1)  records  of  the 
center  element  of  NORSAR  subarray  01A,  used  in  this  study,  were  not 
clipped  and  had  satisfactory  signal  to  noise  ratios;  (2)  accurate  shot 
point  to  surface  velocity  (called  medium  velocity,  hereafter)  values, 
obtained  by  vibrosels  or  similar  techniques,  were  available;  events 
before  1970  could  not  therefore  be  used  (Blandford  et  al.,  1977);  and 
(3)  all  events  were  below  the  water  table  and  confined  to  Yucca  Flats 
and  Pahute  Mesa  regions  of  the  NTS. 

Figure  1  shows  typical  waveforms  for  the  explosions  STRATT  at  Yucca 
Flats  and  POOL  at  Pahute  Mesa.  The  yields  of  the  two  explosions  are 
nearly  equal.  Both  records  have  nearly  the  same  peak  amplitude  values 
but  the  frequency  contents  of  the  Initial  P  wavetraln  appear  to  he  quite 
different.  The  P  codas  on  the  two  records,  however,  do  not  seem  to  be 
significantly  different.  The  time  windows  used  for  spectral  analyses  of 
initial  P  and  P  coda,  as  indicated  in  Figure  1,  are  12.8  sec  and  25.6 
sec  long  with  the  P  window  starting  4  sec  before  the  onset  of  direct  P. 
With  digitization  rate  of  20  samples  per  second,  the  two  windows, 
containing  256  and  512  points,  respectively,  were  tapered  with  a  Parzen 
window  and  Fourier  transformed.  The  two  spectra  were  smoothed  over  4 
and  8  points,  respectively.  The  noise  spectra  were  obtained  In  the  same 
manner  by  using  a  time  window  immediately  preceding  the  Initial  P  window 
and  of  length  12.8  sec.  The  amplitude  spectra  of  observed  signal  (S), 
noise  (N)  as  well  as  S-N,  uncorrected  for  Instrumental  response,  for 
initial  P  and  P  coda  from  STRAIT,  are  shown  In  Figures  2a  and  2b, 
respectively.  The  corresponding  amplitude  spectral  ratio  P/P-coda, 
corrected  for  noise  by  subtracting  the  power  In  noise  from  the  power  In 
observed  signal.  Is  shown  in  Figure  2c.  For  certain  frequency  values, 
the  power  in  noise  exceeded  the  power  in  the  observed  signal.  Points 
for  which  S/N  was  less  than  1.4  were  not  plotted  on  the  P/P-coda  plots. 
Similar  spectra  and  spectral  ratio  P/P-coda  for  MAST  are  shown  in 
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average  P  wave  velocity  from  surface  to  shot  point 


energy  In  the  Initial  P  for  the  two  explosions  with  nearly 
equal  yields  Is  considerably  different  whereas  the  P  coda 
levels  do  not  appear  to  be  much  different. 
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Figure  2a  Results  from  STRAIT  recorded  at  NORSAR.  Displacement  spectra, 
not  corrected  for  Instrumental  response,  of  the  observed 
Initial  P  signal  (S)  denoted  by  the  symbol  +,  noise  (N),  symbol 
0,  for  a  window  of  equal  length,  and  S-N,  symbol  *. 
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Figure  3a  Reaulta  from  MAST  recorded  at  NORSAR,  similar  to  those  for 

STRAIT  in  Figure  2s.  Spectral  nulls  in  initial  P  at  about  l 
and  2  He  are  distinct  and  can  also  be  seen  in  the  spectral 
ratio  plot  of  Figure  3c. 


VARIATION  OF  SPECTRAL  RATTO  P/P-CODA  WITH  EXPLOSION  MEDIUM  VELOCITY 


A  comparison  of  the  spectra  and  spectral  ratios  In  Figures  2  and  3 
shows  seversl  Interesting  features.  P  coda  spectra  appear  considerably 
smoother  then  the  lnltlsl  P  spectra.  The  spectral  holes  observed  In  the 
P  spectra  due  to  cancellation  by  pP  seem  to  have  been  filled  In  the  P 
coda  spectra.  In  fact,  the  spectral  nulls  for  MAST  have  been  estimated 
to  be  at  about  1.0  Hz  and  2.0  Hz  (Shumway  and  Blandford,  1978).  These 
are  clearly  seen  on  the  spectra  of  Figure  la  but  nearly  absent  In  Figure 
3b.  These  spectral  nulls  can  also  be  seen  on  the  spectral  ratio  In 
Figure  3c  as  spectral  minima.  For  STRAIT,  the  spectral  nulls  (see 
Figure  2c)  are  not  as  well  defined  as  for  MAST.  The  spectral  ratio 
P/P-coda,  on  the  average,  decreases  with  frequency  for  STRATT  and 
Increases  with  frequency  for  MAST.  In  other  words,  Initial  P  for 
STRAIT,  an  explosion  In  a  low  velocity  medium  (see  Table  1)  Is  richer  In 
low  frequencies  compared  to  P  coda  whereas  the  opposite  Is  true  for 
MAST,  an  exploalon  In  much  higher  velocity  medium.  The  frequency 
dependence  of  spectral  ratio  P/P-coda  appears  to  be  related  to  the 
medium  velocity  and  probably  other  properties  of  the  near-fleld  medium. 

The  spectral  ratios  P/P-coda  were  also  obtained  for  12  other 
explosions  for  which  the  explosion  medium  velocities  were  known.  The 
least-squares  mean  slopes  on  P/P-coda  plots  (suck  as  those  In  Figures  2c 
and  3c),  are  plotted  against  shot  medium  velocity  In  Figure  4  for  all  14 
explosions.  There  Is  considerable  scatter  among  the  data  points  but  an 
approximately  linear  trend  with  correlation  coefficient  of  0.803  Is 
obtained.  The  large  variation  In  the  ratios  of  P/P-coda  for  the  14  NTS 
explosions  has  to  be  mainly  due  to  differences  in  the  near-source  region 
because  the  receiver  is  same  and  also  because,  for  measurements  at 
teleselsmlc  distances,  the  effects  of  near-receiver  structure  would  tend 
to  be  constant  when  the  ratio  of  P  to  P-coda  Is  tsken. 

Greenfield  (1971)  suggested  that  P  coda  from  explosions  Is 
principally  due  to  near-source  scattering  of  short-period  Rayleigh  waves 
Into  P  waves.  The  large  amplitude  explosion-generated  Rayleigh  waves 
slowly  propagate  outwards  with  appreciable  energy  for  several  tens  of 
seconds  and  a  fraction  of  this  energy  is  gradually  converted  to  P  waves 
and  radiated  downwards  to  teleselsmlc  distances. 


A  possible  qualitative  explanation  for  the  dependence  of  the 
apectral  retlo  P/P-coda  on  the  medium  velocity  can  be  offered  by 
comparing  shots  in  low  and  high  velocity  media.  The  Initial  P,  as 
defined  in  Figure  1,  consists  mainly  of  direct  P,  pP  and  P  waves  derived 
from  conversion  of  near-field  Rayleigh  waves  (called  R  *•  P  waves, 
hereafter).  The  observed  stability  of  P  coda  and  the  variability  of 
first  P  suggest  the  existence  of  small-scale  heterogeneities  In  the 
earth's  Interior  (Alii,  1982).  In  particular,  from  coherency 
measurements  of  near-field  ground  motion  from  an  NTS  explosion, 
McLaughlin  et  al.  (1983)  concluded  that  the  upper  1  km  or  more  of  Pahute 
Mesa  is  heterogeneous  at  scale  lengths  less  than  1  km.  Scattering 
should  therefore  be  greater  for  smaller  wavelength  or  Increase  with 
frequency  and  decrease  with  the  medium  P-wave  velocity.  The  region 
surrounding  an  explosion  in  a  low  velocity  medium  will  act  as  a  high 
scattering  (low  scattering  Q)  region  not  only  because  of  the  lower 
medium  velocity  but  also  because  of  generally  greater  Impedance 
contrasts  which  Increase  the  scattering  efficiency.  Furthermore, 
seismic  energy  from  shots  in  low  velocity  media  may  enounter  low 
inelastic  Q  near  the  source  due  to  lower  hydrostatic  pressure  (Stewart 
et  al.,  1983).  All  these  factors  make  the  initial  P  (comprising  P,  pP 
and  R  -*■  P  waves)  from  an  explosion  in  low  velocity  medium  deficient  In 
high  frequencies.  Due  to  greater  scattering,  the  high  frequency  energy 
is  scattered  away  or  removed  from  the  initial  P  and  some  of  it  may 
appear  later  as  P  coda.  This  is  similar  to  the  apparent  attenuation  of 
a  scattering  medium  described  by  Richards  and  Menke  (1983).  The  Initial 
P  is  therefore  deficient  In  high  frequencies  as  compared  to  P-coda. 

Note  that  the  P  coda  used  in  this  study,  starting  8.8  sec  after  the 
first  P  and  lasting  for  25.6  sec  (see  Figure  1),  is  perhaps  sufficiently 
delsyed  with  respect  to  the  first  P  so  that  its  spectrum  represents  an 
average  over  such  a  large  region  around  the  shot  point  that  the 
scattering  effects  of  localized  low  or  high  velocity  overburden  lose 
their  significance.  The  spectral  ratio  P/P-coda  will  therefore  have  a 
negative  slope  for  explosions  In  low  velocity  medium.  There  will  be 
much  less  scattering  for  explosions  in  high  velocity  medium  so  that  the 
spectral  ratio  P/P-coda  will  have  a  smaller  negative  or  even  positive 
slope.  This  explanation  is  of  course  simplistic  and  purely  qualitative 
because  scattering  is  a  complicated  function  not  onlv  of  the  medium 
velocity  but  also  of  numerous  other  parameters  such  as  the  actual  size, 


shape  and  location  of  the  heterogeneities.  Source  depth  and  the 
presence  of  topographic  features  In  the  source  region  are  also  expected 
to  play  Important  roles.  The  former  strongly  Influences  the  generation 
of  short-period  Rayleigh  waves  (Hudson  and  Douglas,  1975)  whereas  the 
latter  may  strongly  control  the  conversion  of  surface  waves  to 
teleselsmlc  P  (Hudson  et  al.,  1973).  A  detailed  Investigation  using  th 
finite  difference  technique  on  more  realistic  models  of  explosions  In 
various  media  Is  planned  for  the  future. 


P  AND  P  CODA  MAGNITUDES  BASED  ON  SPECTRAL  INTEGRATION 

Using  P  snd  P  coda  windows  specified  earlier,  spectra  corrected  for 
noise  and  lnstruaental  response  were  first  obtained  for  the  1A 
explosions  listed  In  Table  1*  Particle  velocities  integrated  over 
several  frequency  ranges  such  as  0.5  to  3.0  Hr,  0.5  to  2.0  Hr,  and  1.0 
to  3.0  Hr  (referred  to  as  spectral  particle  velocity)  were  obtained  and 
plotted  against  the  known  yield  of  explosions  on  a  log-log  scale.  The 
signal  to  noise  ratio  for  most  explosions  fell  off  rapidly  beyond  3  Hr 
so  that  frequencies  higher  than  3  Hr  were  not  considered  to  be  useful. 

The  frequency  range  O.S  to  3.0  Hr  appeared  to  provide  better  correlation 
with  yield  then  the  other  frequency  ranges.  The  results  for  P  and  P 
coda  are  shown  In  Figures  S  and  6,  respectively.*  The  data  for  P  show 
large  scatter  and  the  Yucca  Plats  and  Pahute  Mesa  events  appesr  to 
separate  out*  For  P  coda,  there  is  no  clear  distinction  between  the 
Yucca  Flats  and  Pahute  Mesa  events.  On  the  basis  of  commonly  accepted 
explosion  source  tine  functlone  (e.g.  von  Seggern  and  Blandford,  1972), 
the  spectrally  Integrated  aagnitude  versus  yield  data  should  show  a  aean 
slope  of  less  than  unity  (see,  for  exaaple,  Blandford,  1976).  An 
atteapt  Is  aade  to  fit  a  slope  of  0.9  to  the  data  in  Figures  5  and  6. 

The  integrated  spectral  aagnttudes  are  subject  to  variations 
slallar  to  those  for  time  domain  measurements.  The  effect  of  pP ,  which 
has  a  modulating  effect  on  the  spectra,  will  in  general  be  somewhat 
smoothed  out  as  compared  to  its  time  domain  counterpart.  P  coda  should 
be  nearly  free  from  the  effects  of  pP.  The  effect  of  attenuation  due  to 
the  upper  mantle  structure  is  also  Important  (Blandford,  1976;  Marshall 
et  si.,  1979)  and  should  equally  apply  to  P  and  P  coda.  A  correction 
for  the  effects  of  attenuation  was  therefore  applied  to  the  spectral 

*  These  two  figures  (as  well  as  Figures  7,8,9  and  10)  do  not  show  the 
actual  data  points  because  the  yield  values  are  classified.  Complete 
figures  are  available  in  Appendix  A  (classified). 


* 
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SPECTRRL  VELOCITY 
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Figure  5 


1.80 

LOG  YIELD 


Spectral  particle  velocity  (arbitrary  units  and  Integrated 
over  the  frequency  range  0.5  to  3.0  Hr)  plotted  against  log 
yield  for  the  initial  P  from  the  fourteen  NTS  explosions. 
Lines  of  slope  of  0.9  are  drawn  through  data  points  for  the 
four  Yucca  Flats  events  (upper  line)  and  the  ten  Pahute  Mesa 
events  (lower  line).  For  actual  data  points,  see  Appendix 


spectral  displacement 


Figure  8  P  coda  Integrated  spectral  displacement  with  t*  and  granite 
aource  function  corrections  (arbitrary  units)  plotted  as  a 
function  of  log  yield.  A  line  of  slope  1.0  through  all  data 
points  has  an  associated  standard  deviation  of  0.06 
magnitude  unit.  For  actual  data  points,  see  Appendix  A. 
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Figure  9 


LOG  YIELD 


P  coda  Integrated  apectral  diaplacenent  with  t*  and  tuff 
aource  function  correctlona  (arbitrary  unita)  plotted  veraus 
log  yield*  A  line  of  alope  1*0  through  13  data  pointa 
(excluding  the  one  with  the  aaalleat  yield)  glvea  a  standard 
deviation  of  0.07  aagnitude  unit.  For  actual  data  points, 
see  Appendix  A. 


SPECTRAL  DISPLACEMENT 

5  4  5  6  789,1  Cf  2  3  H  5  6  789,10* 


LOG  YIELD 


Figure  10  Slallar  to  Figure  9  except  that  the  tuff  source  function 
correction  la  based  on  approximate  yields  derived  from  m^ 
values*  A  line  of  slope  1*0  through  the  thirteen  data  points 
Indicates  a  standard  deviation  of  0.09  magnitude  unit.  For 
actual  data  points,  see  Appendix  A. 


data  by  aaaualng  t  ■  0.4,  the  value  appropriate  for  NTS  ahota  recorded 
at  NORSAR  (Der,  et  al.,  1983).  The  particle  dlaplaceaent  value  for  each 
frequency,  f  was  therefore  aultlplled  by  the  tera  exp  ( n  f  t  ). 

Finally,  a  correction  for  the  aource  function  was  applied  using  the  von 
Seggern  and  Blandford  (1972)  aource  tlae  function  for  the  known  ylelda 
of  the  14  NTS  exploalona.  The  aource  function  for  a  given  yield,  Y  may 
be  expreaaed  aa 


S(w)  -  Y 


[A2<£)2  +  1]1/2 


0)  ?  3/2 

KP  +  1]  ' 


(  «  ) 


(1) 


where  o)  ■*  2  *  f,  Y(«  )  is  a  coupling  tera  Independent  of  frequency,  and 

k  -  k  (5/Y)1^  (2) 

o 

where  A  and  kQ  are  medlua-dependent  constants.  The  source  function 
correction  la  therefore  the  multiplication  of  the  particle  displacement 
values  by  S(o)/S(w).  With  these  two  corrections,  the  Integrated 
spectral  magnitude  versus  yield  plots  should  Ideally  have  a  slope  of 
unity. 


The  14  NTS  explosions  used  In  this  study  were  mostly  In  tuff  and 
rhyolite  (only  one  was  In  alluvium,  see  Table  1).  All  shots  were  below 
the  water  table  so  that,  following  Blandford  (1976),  we  used  the  granite 
reduced  displacement  potential  (RDP)  for  teleselsmlc  P  (l.e.  A  ■  5. OR 
and  kQ  •  16.8).  For  P  coda,  we  first  tried  the  same  granite  RDP.  P 
coda  la,  however,  derived  mainly  from  waves  propagating  upwards  and 
sideways  so  that,  again  following  Blandford's  (1976)  suggestion,  a 
source  function  for  tuff  (A  ■  1.0,  kQ  ■  12.0)  was  also  attempted.  The 
results  for  P  and  P  coda,  with  both  t  and  source  function  corrections, 
are  given  in  Figures  7,  8  and  9. 

Comparing  the  two  sets  of  results,  spectral  integration  of  initial 
P  In  Figure  7  is  similar  to  that  In  Figure  5  but  with  somewhat  smaller 
source  bias  between  the  Yucca  Flats  and  Pahute  Mesa  events.  Considering 
all  14  explosions,  the  standard  deviation,  o  ,  associated  with  fitting  a 
slope  of  1  is  0.14  magnitude  unit  (m.u.).  A  comparison  of  the  results 
for  P  coda  In  Figures  8  and  6  appears  to  show  some  improvement  when  t 
and  source  function  corrections  are  applied.  Moreover,  a  slope  of  unity 
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■eeaa  to  b«  valid  for  the  data  In  Figures  8  and  7.  Obtaining  a  elope  of 

unity  haa  practical  algnlflcance  aince  it  auggeata  that  correctlona  for 

* 

aource  function  and  t  can  appropriately  be  applied  to  P  coda  ao  that 

* 

reaulta  from  several  stations  (each  with  its  own  t  value)  can  be 
properly  coablned  to  further  improve  the  yield  determinations.  The 
standard  deviation  for  the  P  coda  data  in  Figure  R  is  only  0.06  m.u. 

The  corresponding  P  coda  results  based  on  the  source  function  for  tuff 
(Figure  9)  indicate  o  of  0.07  m.u.,  if  the  smallest  yield  explosion  la 
excluded.  The  latter  seems  to  be  an  outlier  perhaps  because  of  reasons 
such  aa  Its  being  the  only  explosion  in  alluvium,  only  34m  below  the 
water  table  and  its  record  at  NORSAR  having  low  signal  to  noise  ratio. 

In  any  case,  the  results  for  P  coda  (Figures  8  and  9)  have  significantly 
less  scatter  than  initial  P,  with  a  ^  0.07  m.u.  Note  that  this  low  value 
of  a  is  nearly  the  same  as  for  yield  estimates  from  regional  Lg  waves 
for  NTS  explosions  when  results  from  several  stations  were  combined 
(Nuttli,  1983). 

In  a  real  situation,  one  would  like  to  use  the  P  coda  from  an 
explosion  to  estimate  its  yield  ao  that  the  value  of  yield  necessary  for 
applying  the  aource  function  correction  will  not  be  available.  We  there¬ 
fore  needed  to  teat  how  good  the  P  coda  estimates  can  be  if  only 
approximate  yields  such  as  those  from  known  m^  values  are  known.  For  NTS 
tuff /rhyolite  events,  Bache  (1982)  suggested  the  following  relationship 
between  magnitude ,  (see  Marshall  et  al.,  1979)  and  yield,  Y 

- w -  -  3.71  +  0.89  log  Y  (3) 

with  o  •  0.09.  One  may  assume  m^  to  be  the  same  as  m^(ISC).  Knowing 
the  m.  value  for  the  13  tuff /rhyolite  explosions  in  Table  1,  we  obtained 

D 

their  approximate  yields  from  equation  (3)  and  applied  the  corresponding 
aource  function  correction  based  on  the  tuff  RDP .  The  results  are  shown 
in  Figure  10  in  which  the  Integrated  spectral  displacement  is  plotted 
against  the  log  of  actual  yield.  Assuming  a  slope  of  unity,  the 
standard  deviation,  a  is  0.09  m.u.,  only  slightly  greater  than  the 
earlier  value  of  0.07  m.u.  Note  that  several  yield  values  derived  from 
the  use  of  equation  (3)  differ  from  the  actual  (classified)  yields  by 
nearly  a  factor  of  2.  Comparing  Figures  9  and  10,  it  seems  that  source 
function  correction  based  on  approximate  yields  works  almost  as  well  as 
when  based  on  actual  yields. 


MAGNITUDES  BASED  ON  SPECTRAL  INTEGRATION  FOR  SOVIET  EXPLOSIONS 


Twenty  three  explosions  from  the  Shagan  River  region  of  the  Eastern 
Kazakh  test  area  of  the  USSR,  with  satisfactory  recordings  at  NORSAR 
were  selected  (Table  2).  P  and  P  coda  magnitudes  were  obtained  using 
the  same  procedure  aa  for  the  NTS  explosions  but  the  frequency  range  was 
Increased  to  O.S  to  5.0  Hz  since  nearly  all  records  showed  good  signal 
to  noise  ratios  for  frequencies  up  to  at  least  5  Hz.  The  correction  for 

attenuation  in  the  upper  mantle  for  Shagan  explosions  recorded  at  NORSAR 

* 

requires  t  ■  0.15  (Der  et  al.,  1983).  Since  the  yield  of  the  Soviet 
explosions  is  not  known,  the  source  function  correction  could  not  be 
applied  in  the  same  maner  as  for  the  NTS  explosions.  An  alternative 
approach,  deacrlbed  below,  was  therefore  used  to  obtain  approximate 
valuea  of  yield  which  were  then  used  to  apply  the  source  function 
corrections. 

Using  an  Improved  general  linear  model  technique  on  a  large  amount 
of  short-period  data,  Blandford  et  al.  (19R3)  deduced  that  a  yield  of 
150  kton  for  a  Shagan  explosion  should  correspond  to  m^  of  6.17. 

Marshall  et  al.  (1979)  defined  a  new  magnitude,  m^  which  correlates  well 
with  yield  for  explosions  at  various  test  sites.  The  explosions  at 
Shagan  are  believed  to  be  In  hard  rock  similar  to  granite  so  that  their 
empirical  relationship  for  explosions  In  salt  and  granite 

niq  -  4.26  +  1  .00  log  Y  (4) 

(see  Figure  6b,  Marshall  et  al.,  1979)  should  be  applicable.  One  may 
also  write 

mQ  -  mb(TS0)  +  C  (5) 

where  C  Is  a  constant.  From  equations  (4)  and  (5) 

log  Y  -  m.(TSC)  +  C  -  4.26  (6) 

D 

Using  the  Information  that  for  Y  -  150  kton,  mb(TRC)  -  6.17,  equation 
(6)  gives  C  ■  0.27  so  that 

log  Y  -  m. (ISC)  -  3.99  (7) 

D 

Knowing  m^dSC)  for  the  Shagan  explosions  (Table  2),  the  yield  can  thus 
be  estimated  and  the  source  function  correction  can  be  applied.  A 
granite  RDP  was  assumed.  Results  for  P  and  P  coda  are  shown  In  Figures 
11  and  12,  respectively,  where  m.(ISC)  Is  plotted  against  the  spectral 

D 

displacement.  There  are  only  14  data  points  for  P  since  P  was  clipped 

for  9  explosions;  P  cods  was  however  available  for  all  events.  The 

results  for  initial  P  (figure  11)  show  a  mean  slope  of  0.75  +  0.05 
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TABLE  2 


SHAGAN  RIVER  EXPLOSIONS  USED  IN  STUDY 


Yield  from 


Date 

Location  (ISC) 

»b(ISC) 

Nuttli 

02 

Nov 

1972 

49.91N, 

78.85E 

6.1 

170 

23 

Jul 

1973 

49.94N, 

78.85E 

6.1 

229 

31 

May 

1974 

49.91N, 

78.91E 

5.9 

62 

16 

Oct 

1974 

49.99N, 

78.96E 

5.5 

20 

27 

Dec 

1974 

49.91N, 

79.05E 

5.6 

63 

27 

Apr 

1975 

49.94N, 

79.02E 

5.6 

35 

29 

Oct 

1975 

49.92N, 

78.91E 

5.8 

33 

25 

Dec 

1975 

50.02N, 

78.86E 

5.7 

93 

21 

Apr 

1976 

49.89N, 

78.83E 

5.3 

17 

09 

Jun 

1976 

49.98N, 

79.07E 

5.3 

21 

04 

Jul 

1976 

49.85N, 

78.97E 

5.8 

112 

28 

Aug 

1976 

49.95N, 

78.98E 

5.8 

49 

29 

May 

1977 

49.86N, 

78.84E 

5.8 

47 

29 

Jun 

1977 

49.96N, 

78.91E 

5.3 

16 

05 

Sep 

1977 

50.05N, 

78.93E 

5.8 

39 

11 

Jun 

1978 

49.88N, 

78.81E 

5.9 

74 

05 

Jul 

1978 

49.84N, 

78.91E 

5.8 

59 

15 

Sep 

1978 

49.91N, 

78 .948 

6.0 

105 

04 

Nov 

1978 

50.03N, 

78.98E 

5.6 

46 

23 

Jun 

1979 

49.89N, 

78.92E 

6.2 

120 

07 

Jul 

1979 

50.05N, 

79 .068 

5.8 

105 

04 

Aug 

1979 

49.86N, 

78.94E 

6.1 

157 

14 

Sep 

1980 

49.94N, 

78 .868 

6.2 

200 

SPECTRAL  DISPLACEMENT 


Initial  P  Integrated  spectral  displaceaent  with  t*  and 
granite  source  function  corrections  (arbitrary  units) 
plotted  as  a  function  of  m.  for  14  Soviet  explosions.  Mean 
slope  is  0.75  +  0.05  and  tne  correlation  coefficient  is 


(i.e.  standard  deviation  of  0.05)  and  correlation  coefficient  of  0.97. 
The  atandard  deviation  of  magnitude  estimate,  <T  is  0.07.  For  P  coda 
(Figure  12),  aean  slope  is  0.70  +  0.03,  correlations  coefficient  ia  0.98 
and  •  0.06.  As  one  would  expect,  correlation  with  a^  is  better  for  P 
coda  than  for  P.  There  is  of  course  no  way  of  checking  how  well  the 
spectral  displacement  values  correlate  with  actual  yields.  It  is, 
however.  Interesting  to  note  that  the  observed  slope  of  about  0.7  for 
both  P  and  P  coda  agrees  well  with  the  slope  derived  for  magnitude-yield 
curves  on  the  basis  of  the  source  function  for  granite  and  t*  ”  0.15 
(see  Figure  7,  Blandford,  1976). 

Using  short-period  Lg  waves  recorded  at  regional  distances  at  a 

number  of  stations,  Nuttli  (1983)  estimated  the  yield  of  all  the  Soviet 

explosions  used  in  this  study;  these  values  are  also  listed  in  Table  2. 

A  plot  of  these  yield  values  versus  coda  spectral  displacement  values 

(Figure  13)  shows  rather  poor  correlation  (correlation  coefficient  of 

0.91).  In  other  words,  our  coda  displacement  values  do  not  correlate 

well  with  Nuttli ' s  yield  values.  It  is,  however,  found  that  the  AFTAC 

classified  m.  values  correlate  better  with  our  coda  displacement  values 
b 

(correlation  coefficient  -  0.97)  than  with  Nuttli’s  yield  values 
(correlation  coefficient  -  0.94).  This  perhaps  means  that  yields  based 
on  P  coda  at  a  single  station  are  at  least  as  accurate  as  those  derived 
from  Lg  recorded  at  a  number  of  stations. 


DISCUSSION  AND  CONCLUSION 


Our  results  from  spectrsl  analyses  of  P  and  P  coda  of  NTS 
explosions  offer  the  exciting  possibility  of  obtaining  source  medium 
Information  from  short-period  teleselsmic  data.  The  positive 
correlation  between  the  spectral  ratio  P/P-coda  and  explosion  medium 
velocity  is  encouraging.  Other  factors  such  as  depth  of  burial, 
subsurface  velocity  and  density  distribution,  topography  are  also  bound 
to  Influence  the  spectra  of  P  and  its  coda  so  that  more  sophisticated 
analyses  incorporating  these  factors  should  lead  to  considerably  more 
information  about  the  explosion  source  and  its  environment. 

Magnitudes  based  on  spectral  integration  of  P  coda  from  NTS 
explosions  provide  a  much  more  reliable  measure  of  yield  than  the 
initial  P.  Furthermore,  the  spectral  displacement  of  P  coda  seems  to  be 
linearly  related  to  yield  after  corrections  for  the  effects  of 
attenuation  in  the  upper  mantle  and  source  function  based  on  knowledge 
of  approximate  yield  (or  mfc)  have  been  applied.  Slope  of  1  on  spectral 
displacement  versus  yield  plots  (on  a  log-log  scale)  suggests  that 
results  from  several  stations  can  be  combined  to  further  improve  the 
yield  determinations.  It  is  also  hoped  that  information  regarding  the 
source  medium,  such  as  that  derived  from  the  spectral  ratio  P/P-coda,  can 
also  be  used  to  further  improve  the  yield  estimates. 

Spectral  Integration  of  P  and  P  coda  from  Soviet  explosions  within 
the  Shagan  River  region  show  m  to  correlate  better  with  spectral 

D 

displacements  from  P  coda  than  those  from  the  initial  P.  There  is  of 
course  no  way  of  verifying  how  well  the  coda  values  correlate  with 
yield.  It  seems,  however,  that  P  coda  from  a  single  station  can 
provide  at  least  as  good  estimates  of  yield  as  those  from  Lg  based  on  a 
network  of  stations. 
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